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ABSTRACT 

• Slowly pulsating B (SPB) stars are g mode pulsators in main-sequence stages with 

I mass ranging from M ~ 3Mq to ^ 8Mq. In this paper, we examine pulsational sta- 

■ bility of low m r modes in SPB stars by calculating fully nonadiabatic oscillations of 

uniformly rotating stars, where m is an integer representing the azimuthal wave num- 
ber around the rotation axis. R modes are rotationally induced, non-axisymmetric, 
oscillation modes, whose oscillation frequency strongly depends on the rotation fre- 
quency f2 of the star. They are conveniently classified by using two integer indices 
I m and I' > \m\ that define the asymptotic oscillation frequency 2mft / [I' {V + 1)] in 

CS| ■ the limit of O ^ 0. We find low m, high radial order, odd r modes with I' — m in 

0^ ■ SPB stars are excited by the same iron opacity bump mechanism that excites low 

frequency g modes of the variables, when the rotation frequency f2 is sufficiently high. 
No even r modes with low m are found to be pulsationally unstable. Since the surface 
pattern of the temperature perturbation of odd modes is antisymmetric about the 
equator of the star, observed photometric amplitudes caused by the unstable odd r 
modes with V = m are strongly dependent on the inclination angle between the axis 
O \ of rotation and the line of sight. Applying the wave-meanflow interaction formalism 

to nonadiabatic r modes in rapidly rotating SPB models, we find that because of the 
r(f> component of the Reynolds stress and the radial transport of the eddy fluctuation 
of density in the rotating star, the surface rotation is accelerated by the forcing due to 
the low V = m unstable r modes. We suggest that the amount of angular momentum 
. redistribution in the surface region of the stars can be comparable to that needed to 

^ \ sustain decretion discs found in Be systems. 

We make a brief comparison between nonadiabatic r mode calculations done with 
and without the traditional approximation. In the traditional approximation, the local 
horizontal component of the rotation vector is ignored in the momentum conser- 
vation equation, which makes it possible to represent the angular dependence of an 
oscillation mode using a single Hough function. We find that the oscillation frequencies 
of low m r modes computed with and without the traditional approximation are qual- 
itatively in good agreement. We also find that the pulsational instability of r modes in 
the traditional approximation appears weaker than that without the approximation. 
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1 INTRODUCTION 

Slowly pulsating B stars (hereafter SPB stars) are g mode pulsators (e.g., Waelkens 1991, Waelkens et al 1998). Low frequency 
g mode oscillations in SPB stars are excited by the iron opacity bump mechanism, which also excites high frequency p modes 
in 13 Cephei stars (Dziembowski & Pamyatnykh 1993; Dziembowski, Moskalik, & Pamyatnykh 1993; Gautschy & Saio 1993, 
see also Kiriakidis, ElEid, & Glatzel 1992). Savonije (2005) and Townsend (2005b) have recently shown that low m r modes in 
SPB stars are also excited by the same opacity bump mechanism that excites g modes, if rapid rotation is assumed. Although 
SPB stars are not necessarily rapid rotators, a fraction of SPB stars are rapidly rotating (e.g., Aerts et al 1999). Rapidly 
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rotating SPB stars will be good observational targets for both g mode and r mode oscillations, which can be used as a probe 
to investigate the interior structure of the variables. 

R modes arc rotationally induced, non-axisymmctric, low frequency oscillation modes whose toroidal component of the 
displacement vector is dominant over the horizontal and the radial components, and whose oscillation frequency is proportional 
to the rotation frequency f2 (e.g., Papaloizou & Pringle 1978; Provost, Berthomieu, & Rocca 1981; Saio 1982; Berthomieu 
& Provost 1983). R modes arc retrograde modes in the corotating frame of the star, but are progradc when soon in an 
inertial frame. In the limit of O ^ 0, the displacement vector is well represented by a single spherical harmonic function 
YIT' and the oscillation frequency oj observed in the corotating frame of the star tends to the asymptotic frequency given by 
u}c{m,l') = 2mCl/[l' {I' + 1)] with I' > \m\. R modes are therefore conveniently classified by using the two integer indices m and 
/'(> |m|) that define the asymptotic freqyency Lodm, I'). For a neutrally-stratified (isentropic) star, in which the Brunt- Vaisala 
frequency N vanishes identically in the interior, the odd r mode with I' = m that has no radial nodes of the eigenfunction is 
the only r mode we can find for given m and (sec, c.g, Saio 1982, Lockitch & Friedman 1999, Yoshida & Lee 2000a). The 
restoring force for the r mode is exclusively attributable to the Coriolis force. For a stably-stratified star, however, buoyant 
force comes into play as an additional restoring force for the low frequency oscillations. For any combination of m, l'{> \m\), 
and fl, there exist an infinite sequence of r modes that differ in the number of radial nodes of the eigenfunction (e.g., Yoshida 
& Lee 2000b). In this paper, we shall call r modes in stably-stratified stars buoyant r modes. 

For r modes to be pulsationally excited in SPB stars, rapid rotation of the star is required (Savonije 2005, Townsend 
2005b). Expecting a fraction of SPB stars arc rapidly rotating, and noting that some Be stars are found in the SPB domain in 
the HR diagram, it may be reasonable to expect that there exists a link between a class of Be stars and rapidly rotating SPB 
stars. Be stars are characterized by emission lines, which are thought to have their origin in the gaseous discs surrouonding 
the central B stars. Time variability of the emission lines in strength and shape as well as their appearance and dissapearance 
are believed to be closely related to the dynamical status of the gaseous discs aurrounding the stars (see a recent review by 
Porter & Rivinius 2003). Although discs around Be stars are indispensable for observed Be phenomena, no definite physical 
mechanisms responsible for disc formation have not so far boon indentified. Considering that both g modes and r modes 
become pulsationally unstable in rapidly rotating SPB stars, it is tempting to speculate these low frequency oscillation modes 
play a role in angular momentum transfer in the interior that leads to disc formation around the stars. To investigate angular 
momentum redistribution in rotating stars, in a series of papers Ando (1981, 1982, 1983, 1986) has applied to stellar pulsations 
the wave-meanflow intercation formalism, which has been extensively used in geophysics and fluid mechanics (e.g., Lighthill 
1978; Craik 1985; Pedlosky 1987; Andrews, Holton, & Leovy 1987). Carring out a numerical simulation, Ando (1986) has 
argued that quasi-periodic vacillation phenomena in Be stars occur as a result of the interaction between the rotation and two 
prograde and retrograde global oscillation modes. Quite similar calculation has also been done to simulate the quasi-biennial 
oscillation in the earth's atmosphere (Plumb 1977; see also Lindzen & Holton 1968). We think it worthwhile to apply the wave- 
meanflow interaction formalism to pulsationally unstable low frequency modes in rapidly rotating SPB stars, speculating that 
unstable prograde and retrograde g modes are responsible for quasi-periodioc vacillations, and unstable r modes are working 
as a stable angular momentum supplier to the surrounding discs. 

In the analyses of nonadiabatic oscillations of SPB stars, Savonije (2005) and Townsend (2005a, b) exclusively used the 
traditional approximation, neglecting the local horizontal component of the rotation vector fl in the momentum conservation 
equation (e.g., Eckart 1960). In the traditional approximation, separation of variables between spherical polar coordinates 
(r,6,(j>) becomes possible and the governing equations for oscillations in rotating stars are largely simplified (e.g., Lee & 
Saio 1987a). The traditional approximation has been extensively employed in geophysics to study waves propagating in an 
atmospheric or ocean fluid layer that is very thin campared to the radius of the earth. The approximation is thought to 
be valid for low frequency oscillations with uj « N so that the horizontal component of the fluid velocity is much larger 
than the radial component (e.g., Eckart 1960). When we apply the traditional approximation to global oscillations of rotating 
stars, however, validity of the approximation is not very clear, since the condition uj « N may not always be satisfied in the 
entire interior, for example. In this paper, we calculate nonadiabatic low m buoyant r modes of SPB stars without applying 
the traditional approximation, in order to examine their pulsational stability. We also apply the wave-meanfiow interaction 
formalism to nonadiabatic r modes to see if acceleration of the surface rotation is possible as a result of the forcing. The result 
of the stability analysis is given in §3 and the wave-meanflow interaction for the r modes is discussed in §4. A brief comparison 
between nonadiabatic r mode calculations with and without the traditional approximtion is given in §5. We briefly describe 
our numerical method in §2, and conclusions are summarized in §6. 



2 NUMERICAL METHOD 

The numerical method of calculating nonadiabatic oscillations of uniformly rotating stars is the same as that given in Lee 
& Saio (1987T3, 1993). We employ a linear theory of oscillations to describe small ampUtude pulsations of rotating stars. 
We assume that the background state of rotating stars is axisymmetric about the rotation axis, and that the time t and 
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azimuthal angle 4> dependence of perturbations to the background state is given by exp[i{at + m4>)], where a is the oscillation 
frequency in an inertial frame, and m is an integer representing the azimuthal wavenumber around the rotation axis. Under 
these assumptions, linearized hydrodynamic equations are given by 



p' + v-(pO = o, 

\sjjpT5s = (pe-V-F)', 
and 

— = 1 1 hCCTl 1 — , 



-Vp' - p'V* - pV$', 



(1) 
(2) 
(3) 
(4) 

(5) 



P P Cp 

where p is the mass density, p is the pressure, T is the temperature, s is the specific entropy, Cp is the specific heat at constant 
pressure, $ is the gravitational potential, F is the energy fiux vector, e is the nuclear energy generation rate per unit mass, 
f2 is the vector of angular velocity of rotation, 



u! = a + rn^l 

is the oscillation frequency observed in the corotating frame of the star, and 
dlnp\ _ fdlnp" 



ri = 



\d\npj^' ^ \dlnTjp- 



(6) 
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Here, ^ is the displacement vector, and the quantities attached by a prime (') and S represents their Eulerian and Lagrangian 

perturbations, respectively. The total energy flux is given as f = i^conv + J'rad, whore i^conv is the convoctivo energy flux, 
which is calculated using a mixing length theory of turbulent convection, and i^rad denotes the radiative energy flux given by 
ac 
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where k is the opacity, and a and c denote the radiative constant and the velocity of light, respectively. For the pertubation 
of the convective flux, we assume 6{V ■ Fconv) = for simplicity. 

We apply boundary conditions both at the surface and at the center of the star. The surface boundary conditions we 
use at r = ii are 6p = 0, S{Lr — 47rr^asBT*) = 0, and the condition that is continuous at the stellar surface, where 
Lr = Awr'^Fr is the luminosity at r, and usb is the Stefan-Boltzmann constant. The inner boundary conditions at r = 0, 
on the other hand, come from the regularity condition of the mechanical variables at the center and from the condition that 
oscillations are nearly adiabatic, Ss/cp ~ 0, in the deep interior. 

In this paper, we employ a spherical polar coordinate system (r, 9, 0) whose origin is at the center of the star and the axis 
of = is along the rotation axis. For a given m, the displacement vector ^(r, 6, 4>, t) is given by a truncated series expansion 
in terms of spherical harmonic functions Vi'"(9, (j)) with different Z's as 



3max 



Si^{r)Y{;{e,<j>)e'' 



6 ^ 

rp / J 



Hi.{r) d 



sin 6 d4> 



and the pressure perturbation, p'(r,6,4>,t), for example, is given as 



(9) 



(10) 



(11) 



(12) 



where Ij = |m| + 2(j — 1) and I'j = Ij + 1 for even modes and Ij = |m| + 2j — 1 and I'j = Ij — 1 for odd modes where 
j = 1, 2, ■ ■ ■ . j'max- Note that the function p' is symmetric (antisymmetric) for oven (odd) modes about the equator of the 
star. Substituting the expansions (9) to (12) into linearized hydrodynamic equations (1) to (5), we obtain a set of coupled linear 
ordinary differential equations of a finite dimension jmax for the radial expansion functions. With the boundaxy conditions at 
the center and at the surface of the star, we solve the set of linear differential equations of a finite dimension as an eigenvalue 
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problem for the oscillation frequency a; (Lee & Saio 1987b, 1993). Oscillation modes with wi = Im(ci;) < are pulsationally 
unstable while those with wi > stable. Note that we do not include terms associated with the centrifugal force, that is, 
setting / = in the oscillation equations given in Loo & Saio (1987b, 1993). 

If a; and ^ are respectively the eigenfrequency and the eigenfunction of a mode governed by the oscillation equations (1) 
to (5), multiplying equation (1) by which is the complex conjugate of the displacement vector 4, and integrating over the 
whole volume of the star, we obtain 

J^E + 2ujF-G = Q, (13) 
where 

E = j pC-idV, (14) 

' pe-{nxi)dv, (15) 



G = j C ■ (Vp + p'V* + pV*') dV. 



(16) 



Note that E and F are real quantities, because (^* • =£,*■£, and [i^* • (fl x ^)]* = i^* • (f2 x ^). Since the term G can be 

rewritten as 



G= dV 



ri p Vp p J P 47rG Cp 



+ JdS- (^Cp + Pr$' + 1 , (17) 



substituting u) = + iwi into equation (13), we obtain 



where we have assumed that Vp is parallel to Vp and that the surface terms can be ignored because of the outer boundary 
conditions employed in this paper. For oscillations of rotating stars, we may define a work integral as 

w(r) = — — i -rr i Im ( — QT(5p^^ ) 47rr^dr, (19) 



2{lo^E + F) 

so that we have w{R) > for unstable modes having toi < 0. We can use w{R) as a check on numerical consistency in 
nonadiabatidc calculations. As f2 — » 0, F — > 0, and the work integral (19) reduces to that for non-rotating stars. We note that 
the region in which dw{r)/dr > is a driving region for the oscillation while the region of dw{r)/dr < a damping region. 



3 MODAL PROPERTY OF R MODES IN SPB STARS 

Main sequence models used for pulsation calculation arc computed with a standard stellar evolution code, into which we 
incorporate the opacity tables computed by Iglesias, Rogers, & Wilson (1992), and Iglesias and Rogers (1996). Figure 1 shows 
evolutionary tracks, from the ZAMS to early hydrogen shell burning stages, of stars having M = SM© to M = SMq for 
X — 0.7 and Z — 0.02, whore the circles arc SPB stars tabulated in Waelkens et al (1998) and the filled circles indicate 
rapidly rotating SPB stars given in Aerts et al (1999). 

We classify buoyant r modes using two integer indices m and l'{> \m\) that define the asymptotic frequency u!c{m,l') = 
2mQ/[l' {I' + 1)] for a given fl. Note that the mixed modes discussed by Townscnd (2005b) correspond to the odd r modes with 
= m in our classification. For given (m, /'), we introduce an integer index n to indicate the radial order of the mode, and we 
write Tn and/or uj„ to denote the n th radial order r mode. The integer n usually corresponds to the number of radial nodes 
of the eigenfunction. In this paper wc exclusively assume m > 0. Then in our convention, the r mode frequency w observed in 
the corotating frame is positive, and the frequency a in an inertial frame is negative. For given (m, /') and fl, the asymptotic 
frequency ujc also sets the upper limit to the r mode frequencies in the sense that ujcim,!') > wo > wi > a;2 > • • • (but see 
also Yoshida & Lee 2000a). We may call ro with loq the fundamental r mode. 7? modes with — (^c have almost solenoidal 
displacement vector, so that V • ^ = —5p/p ~ and hence Sp/p ~ in the interior. 



3.1 Pulsational stability of buoyant r modes 

We calculate nonadiabatic low m r modes of uniformly rotating staxs in main-sequence stages for M = 3Mq to 8Mq for 
X = 0.7 and Z = 0.02. As a typical example of SPB stars, we mainly discuss the case for the 5Mq main-sequence models. 
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The results for odd r modes with T = m = 1 and 2 for = Q./^/GMJW' = 0.4 are given in Figures 2 and 3, respectively, 
where the dimensionless oscillation frequency wr = ojb./ \J GM/R^ is plotted as a function of LogTes with Tgs being the 
effective temperature from the ZAMS to the TAMS. Here, the small thin and large thick circles indicate pulsationally stable 
and unstable modes, respectively. Note that the dimensionless limiting frequencies CL!c{m,l') for the r modes with I' = m = 1 
and 2 are respectively ujc = 0.4 and 0.2666 for Q, = 0.4. Figure 4 is the same as Figure 2 but for Q, = 0.2, for which ujc = 0.2. 
Since the oscillation frequency a; of r modes is essentially proportional to fl, the angular rotation frequency as large as H^.l 
is necessary for the r modes to be effectively excited by the opacity bump mechanism. As shown by these figures, only high 
radial order rn modes having n >> 1 become pulsationally unstable. As the star evolves from the ZAMS to the TAMS, 
the frequency range and the number of unstable r modes of the SAf© star gradually increases and the frequency spectrum 
becomes denser, reflecting the developement of the fj, (mean-molecular-weight) gradient zone just outside the convective core. 
Note that all the r modes are stabilized before the star goes into shell hydrogen burning stages. Compared with the case of 
n — 0.2, the width of Toff for unstable r modes is wider for Cl = 0.4 in the HR diagram, suggesting that rapid rotation favors 
the excitation of r modes. 

As found from Figures 2 to 4, in the course of evolution from ZAMS, the fundamental ro modes with I' = m, which behave 
like an interface mode having the maximum amplitude at the interface between the convective core and the radiative envelope, 
become unstable because of the e mechanism (see also Townsend 2005b). The growth rates r] = —loi/lor of the modes, however, 
are of order of 10~^^ to lO"'^^, and the unstable fundamental ro modes may not have important observational consequences. 

In Figure 5, we plot, as an example, the work integral ■w{r) for a unstable odd r2o mode with /' = m = 1 of a 5Mq 
model with LogTes = 4.188 for fl = 0.4. In this figure, we also plot /tad = (91n/t/91np)ad versus r/R for the model. Strong 
excitation for the rgo mode occurs in an outer envelope region where dnjdr is positive, and overcomes damping effects in 
the deep interior. In Figure 6, the growth rate r\ of unstable r modes with /' = m = 1, 2, and 3 of the hMq, model with 
LogTeff = 4.188 is plotted versus the inertial frame oscillation frequency |c7-| for H = 0.4 (O.OllSmHz), where the two filled 
wedges on the horizontal axis indicates the rotation frequencies f2 and 20. The r modes with I' = m = \ have on average 
larger growth rates, compared to those of the r modes with I' = m = 2 and 3. 

In Figure 7, we plot of even r modes with = m -t- 1 = 2 versus LogTeff for Cl = 0.4. No pulsationally unstable r 
modes with even parity are found. This may be attributable to the fact that the even r modes with Z' = m -|- 1 = 2 for a given 
n have frequencies by about a factor of 3 smaller than those of the odd r modes with Z' = m = 1. In fact, we have tried to 
find unstable even r modes assuming Q, = 0.8, but failed. Figure 8 depicts the work integral w{r) of an even r2o mode with 
r = m -h 1 = 2 of the 5Mq model with LogTeff = 4.188 for fj = 0.4. The damping effect in the deep interior is by a factor of 
2 to 3 larger than the excitation near the surface. 

It may be instructive to mention a result for low m r modes of more massive stars, since the instability for low frequency 
oscillation modes weakens as the stellar mass increases. Figure 9 shows cDr of odd r modes with Z' = m = 1 of 8Mo main- 
sequence models from the ZAMS to the TAMS for H = 0.6, where the small thin and large thick circles indicate pulsationally 
stable and unstable modes, respectively. For Q, = 0.6, high radial order r modes with I' = m = 2 and 3 become unstable too 
in the course of evolution, although only the odd r modes with Z' = m = 1 become unstable for CI = 0.4. 



3.2 Eigenfunctions of buoyant r modes 

Let us first show some examples of the eigenfunctions (expansion functions) of low m r modes of the 5Mq ZAMS star, which 
consists of the convective core and the radiative envelope that contains two thin convection zones in the outer part of the 
envelope. Figure 10 depicts the eigenfunction a::Re(iT;/ ) of odd r modes with Z' = ?n = 1 (panel a) and of even r modes with 
I' = m + 1 = 2 (panel b) for Q = 0.4, where x = r/R, and the amplitude normalization is given by Re(S';j) = 1 at the surface. 
Note that in the vicinity of the outer boundary of the convective core, the even r modes have a node, which is not necessarily 
apparent in the figure. Since V — Vad is assumed to vanish in the convective core, iT;/^ is almost constant there (see, e.g., Saio 
1982). Although the amplitude of the ro mode has the maximum just outside the convective core, the amplitude maximum 
shifts to the surface, as the radial order n and the number of radial nodes in the radiative envelope increase. It is important 
to note that for the even r modes with Z' = m + 1 = 2, the isentropic convective core is an evanescent region for the oscillation 
and no fundamental mode that has no radial nodes of the eigenfunction Re(iT;/^) exists. We find this is also the case for both 
even and odd r modes if Z' > \m\. 

Figure 11 shows the eigenfunction a;Re(iT;/ ) of odd r modes with Z' = m = 1 of a 5Mq model with LogTog = 4.188 for 
n = 0.4, where the long dash-dotted line, short dash-dotted line, dotted line, and solid line indicate ro, ri, rio and r2o modes, 
respectively, and the amplitude normalization is given by Re(Sji ) = 1 at the surface. This model has the fj, gradient zone just 
outside the convective core. The fundamental ro mode, which is shown in the inlet, behaves like an interface mode and the 
amplitude has the sharp maximum at the interface between the convective core and the /^-gradient zone. We also note that 
the high radial order r modes have very short radial wavelengthes in the /^-gradient zone, which leads to the dense frequency 
spectrum. 
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4 FORCING ON ROTATION BY BUOYANT R MODES 

In the wave-meanflow interaction formalism, every physical quantity of a rotating star is decomposed into the zonally- averaged, 
axisymmetric part and the non-axisymmetric, residual part, and we regard the former as the quantity in the meanflow and 

the latter as the one representing waves (oscillations) (e.g., Pedlosky 1987, Andrews, Holton, & Leovy 1987). The forcing 
equation for stellar rotation may be given by (e.g., Lcc & Saio 1993; see also Ando 1983) 

^j(nt) = -\-^^^'^^ (rsin^Re {pv'^v* +V4,pVr))] - i (Re (imp *!>')> , (20) 
where 

__ _ 1 _ 1 

j{r,t) = {r sine pv^) , = — J v^d(p, p=—j pd<j>, (21) 

and 

{■■■) = l£ dOsinei---), (22) 
and we rewrite the forcing equation (20) as 

To(r) Ti(r) T2(r) T3(r)' ^'^^^ 
where r's are time scales of local variation of the meanflow defined as 
1 dj{r,t) 



To 



Tl 



T2 



j{r,t) dt 

1 1 1 



2 j{r, t) dr 

I l ia 

2 j{r,t) dr 
1 1 



(24) 

r^(r sin 6'Re(pi;^v;*))], (25) 
(rsin^Re (w^p'wr*))], (26) 



2j(r,t)^^'('"'^'**')^' ^'^^ 

and positive to stands for acceleration of the rotation. Note that v' = iuj^ for uniformly rotating stars. The importance of 
non-conservative (e.g., non-adiabatic) effects accompanied with the oscillation will be apparent by rewriting the Eulerian 
formulation into a Lagrangian one (e.g, Ando 1986b; see also Andrews & Mclntyre 1978). The time derivative of the total 

angular momentum, J, may be given by 



= / 4irr^^dr=-2TTR'^(Rsm9Re{pv'^v!; +v^pv!;))- f 27rr^ (Re (imp'*l>')) dr, 
Jo Jo 



(28) 



where the first term on the right hand side stands for the angular momentum luminosity, and the second term is the gravita- 
tional torque. 

In Figure 12, we plot t^^{i~) with j = 0, 1, and 2 versus r/R for a unstable I' = m = 1 r2o mode of the 5Mq model 
with LogTeff — 4.188 for Q = 0.4, where the solid line, dotted line, and dashed lines represent T(j~^(r), Tj~^(r), and T^^(r), 
respectively, and the amplitude normalization for the mode is given by Kc{Si-^{R)) = 1. For self-excited oscillations, the term 
T^^ is negligible compared to the other two terms and is not shown in the figure. Although the term T2^{r) was ignored 
in the analysis by Ando (1983), the term for low frequency oscillation modes can be more important than the term Tj"^ in 
the surface region because of the boundary condition 5p — Q we apply at the surface. Since Tf^^{r) is positive in the region 
of r ^ R, the unstable r20 mode contributes to accelerating the surface rotation. This acceleration of the surface rotation 
is brought about by many unstable I = m r modes in SPB stars as shown by Figure 13, in which Tq^{R) and j/10^* for 
odd r modes with /' = m = 1 of the ^Mq model for Q, = 0.4 are plotted versus the dimensionless oscillation frequency cDr. 
Here, the circles and squares stand for Tq^{R) and j/10^"', respectively, and the filled and open symbols indicate pulsationally 
unstable and stable modes, respectively. Note that to calculate the quantities Tq^{R) and j, we have employed the amplitude 
normalization given by Re(S';j {R)) = 1 for the modes. We find that Tq^{R) is positive and J is negative for the pulsationally 
unstable I' = m = 1 r modes of the model. Negative J means a loss of angular momentum from the star, or equivalently a 
positive angular momentum luminosity through the surface of the star. 

The numerical values of J in Figure 13 can be compared with observationally estimated ones. The amount of angular 
momentum supply in a unit time necessary to sustain a decretion disc around the central star in a Be system may be estimated 
as (e.g., Lee, Saio, & Osaki 1991) 
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where Md denotes the mass decretion rate through the disc, and M and R are the mass and radius of the central star, and 
an estimation of Md ~ W~^°MQ/yr given by Okazaki (2001) has been used. Since the quantity J defined in equation (28) 
is proportional to the sqauro of the oscillation amplitudes, if we consider the typical value of Jd is of order of 10^^ to 10^^ 
for a decretion disc in Be stars, the oscillation amplitudes of order of Re{Si^{R)) = 0.01 ~ 0.1 may be able to produce J 
comparable to the observational estimations Jd ■ 



5 COMPARISON WITH THE TRADITIONAL APPROXIMATION 

It may be appropriate here to make a brief comparison between numerical results obtained with and without the traditional 

approximation for r mode oscillations of rotating stars. In this section, for both calculations with and without the traditional 
approximation, we consistently employ the Cowling approximation, neglecting the Eulerian perturbation of the gravitational 
potential. 

The traditional approximation has been used by several authors (e.g., Ushomirsky & Bildston 1998, Townscnd 2005a,b; 
Savonije 2005) to numerically investigate low frequency oscillations of SPB stars. Under the traditional approximation, we 
ignore the local horizontal component of the rotation vector SI in the momentum conservation equation. Adiabatic oscillation 
equations in the Cowling approximation then becomes sclf-ajoint, and separation of variables becomes possible between 
spherical polar coordinates (r, 9, (f>) so that an oscillation mode of a rotating star can be represented by 

p\r,e,<t>,t) = Vkrr.ni:r)Hkm{cose)e'"'^+"'\ (30) 

where n denotes the radial order, and the Hough function ^/^^(cos^) with integer indices k and m is the solution to Laplace's 
tidal equation and depends on the ratio 2Q,/oj (e.g., Bildsten, Ushomirsky, & Cutler 1996; Lee & Saio 1997). This separation 
of variables greatly simplifies the problem of analyzing numerically and analytically low frequency oscillations of rotating 
stars, although couplings between oscillation modes associated with different Hough functions cannot be treated under the 
approximation. In geophysics, the traditional approximation has long been used to discuss low frequency atmospheric or ocean 
waves that satisfy N/u) » 1 and have the horizontal velocity much larger than the vertical (radial) velocity (e.g., Eckart 
1960). Although atmospheric or ocean waves discussed in geophysics are propagating in a layer much thinner compared to 
the radius of the earth, low frequency modes of a rotating star in general have eigenfunctions extending from the center to 
the surface of the star and do not always satisfy in the interior the conditions needed for the tradiational approximation to 
be valid. 

A typical difficulty in applying the traditional approximation to stellar pulsations appears in spurious singular behavior 
of the eigenfunctions at the center for certain groups of low frequency modes. In the tradiational approximation, the governing 
equations for adiabatic oscillations may be given by (e.g., Lee & Saio 1987a) 

d£i { V n,\ .. . f ^km V ' 
dr 



= [ciu>^ -rA)zi + {l-rA-U)z2, (32) 



dZ2 

dr 

where the variables zi and Z2 respectively correspond to ^r/r and p' /pgr, and Xkm is the separation factor, which depends 
on the ratio !>" s 20, /u! for given integer indices m and k, and 



_ din Mr dlnp ^ dlnp 1 dlnp ^ {r/Rf 

dlnr' dlnr' dlnr Tidlnr' '^^ M,/M' ^ ' 

and Mr = Airr^pdr. Here, we have used the notation for Xkm employed by Lee & Saio (1997). For A: > 0, Xkm > tends to 
(|m| + fc) (|m| + + 1) as z/ — > 0, while for fc < 0, we have branches of u, in which Xkm is negative. R modes are associated with 
Xkm with negative k that is positive for cj < 2mQ/[l'{l' + 1)] where I' = m + \k + 1\. Assuming both zi and Z2 in equations 
(31) and (32) are proportional to r'^ in the vicinity of the stellar centre, we obtain the regularity condition at the centre for 
the functions zi and Z2 given by 

/3 = (-5 + Vl + 4:Xkm)/2 = i-2>-l, (34) 
which is rewritten as 

Xkm = !{! + I) > 2. (35) 

As suggested by Figure 1 of Lee & Saio (1997), prograde modes associated with k = and m = 1 and r modes with fc < — 1 

except for high radial order fc = — 1 r modes do not satisfy the condition (35), and they show spurious singular behavior of 
the eigenfunctions zi and 22 at the center that does not exist without the traditional approximation. 

To avoide this mode-dependent, awkward singular behavior we meet when numerically integrating the oscillation equations 
(31) and (32), Townsend (2005a,b) introduced a new set of dependent variables qi and 92 defined by 



8 U. Lee 



Z2 = q2X' 



(36) 



to obtain 



(37) 



= (cio)" -rA)qi + (l-rA-U- /?) 92, 



(38) 



where x = r/R. Obviously, the variables qi and 92 are always regular at the stellar center, and we have no trouble to solve 
numerically equations (37) and (38) to obtain cigonfroquoncios lu of low frequency modes. 

Under the traditioanl approximation, Townsend (2005a,b) numerically solved a set of coupled linear ordinary differential 
equations derived by using regularized dependent variables like qi and 92 for adiabatic and non-adiabatic oscillations. Here, 
we calculate non-adiabatic r modes of the 5Mq ZAMS model with Z = 0.02, following the formulation given by Townsend 
(2005a), and we compare the results to those computed without the traditional approximation. Figure 14 shows ujn (panel 
a) and the growth rate \r)\ (panel b) of odd r modes with I' = m = 1 versus the radial order n for = 0.4. Prom panel (a), 
we find that the two calculations for cDr arc qualitatively in good agreement with each other, but it is also apparent that the 
difference in ujn between the two calculations for a given radial order n becomes comparable to the difference in wr between 
two consecutive radial order modes as the radial order n increases. Panel (b), on the other hand, shows that the growth 
rates of unstable r modes obtained in the traditional approximation are systematically smaller and the instability band in 
the radial order n is narrower, compared to those obtained without the approximation. We also find that the fundamental ro 
mode is unstable in the traditional approximation, but it is not without the approximation. This difference in the stability of 
the fundamental r modes may be attributable to the singular behavior of the cigcnfunctions in the traditional approximation, 
since the large amplitudes of the eigenfunctions may wrongly strengthen the instability in the convective core. 

Pigure 15 plots (Dr versus the radial order n for even r modes with T = m + 1 = 2 for = 0.4, where the small filled 
and large open circles stand for the results with and without the traditional approximation, respectively. Although the two 
calculations of wr are in reasonable agreement, it is interesting to note that in the traditional approximation there exists the 
fundamental ro mode with even parity, which does not exist if we lift the traditional approximation. 

In Figure 16, we plot Ke{Xkm) versus the radial order n of r modes, where the filled and open circles respectively stand 
for odd r modes with I' = m = 1 and even r modes with = m + 1 = 2 for H = 0.4. Note that the former and the latter 
corresponds to (fe,m) = (—1, 1) and (—2, 1), respectively, since I' = \m\ + |A; + 1| for r modes. The figure clearly shows that 
the low radial order r modes with I' = m = 1 and all the r modes with I' = = 2 Ao not satisfy the regularity condition 

Re(Afem) > 2 for the functions zi and Z2 at the center. As an example. Figure 17a depicts the eigenfunctions xKe(zi) and 
a;Re(yS;i) for the fundamental r modes with I' = m = 1, while Figure 17b shows the eigenfunctions for no modes with 
I' = m = 1, whore f2 = 0.4, and the functions zi and Si^ arc obtained with and without the traditional approximation, and 
the amplitude normalization for z\_ and Si^ is given by Re(2;i) = 1 and Re(5'zi) = 1 at the surface. The function xz\_ of the 
fundamental mode is singular at the stellar center, but xSi^ is not. On the other hand, no singular behavior appears for the 
high radial order rio modes with T = m = 1 oven in the traditional approximation. 

It seems pulsational instability of r modes appears weaker in the traditional approximation. This trend can be seen more 
clearly in pulsational stability of low m r modes of more massive stars as shown by Figure 18, in which wr of T = m = 1 
r modes computed for 8Mq main-sequence models in the traditional approximtioan is plotted versus LogTcff for = 0.6, 
where the samll open circles and large filled circles stand for stable and unstable modes, respectively. Camparing Figure 18 
with Figure 9, we find the instability domain in the traditional approximation is largely shrinked and the blue edge of the 
instability strip has shifted to the lower T^s in the HR diagram. 



6 CONCLUSIONS 

In this paper, we have calculated nonadiabatic r modes of SPB stars with mass ranging from M = SM© to M = 8M0 , and 

found that low m, odd r modes with I' = m become pulsationally unstable if we assume rapid rotation of the stars, confirming 
the numerical results obtained in the traditional approximation by Savonije (2005) and Townsend (2005b). No low m, even 
r modes with I' = m + 1 are found pulsationally unstable, which contradicts Savonije (2005) who found unstable even r 
modes. A possible reason for the contradiction may be that Sovonije (2005) calculated forced oscillations for r modes of SPB 
stars, which are not necessarily equivalent to free r mode oscillations computed in this paper. Photometric amplitudes of 
the variability caused by unstable odd r modes with I' = m will be strongly dependent on the inclination angle between the 
rotation axis and the line of sight, because the surface pattern of the temperature perturbation is antisymmetric about the 
equator of the star. Applying the wave-meanfiow formalism to pulsationally unstable r modes in SPB stars, we have found 
that the unstable r modes have a contribution to accelerating the surface rotation of the star. Note that it is also important 
to examine the mechanism for unstable g modes. 
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We have made a brief comparison between r mode calculations done with and without the traditional approximation. 
We found that the oscillation frequency wr of r modes obtained in the traditional approximation are qualitatively in good 
agreement with those without the traditional approximation. The growth rates of unstable r modes in the traditional ap- 
proximation are systematically smaller than those obtained without the traditional approximation, that is, the pulsational 
instability appears weaker in the traditional approximation. For r modes associated with Xkm < 2, the eigenfunctions obtained 
in the traditional approximation show spurious singularity at the center, which may affect the pulsational stability, particularly 
for the low radial order modes. Note that we need to make a comparison for the case of g modes to have a good understanding 
about the validity of the traditional approximation. For low frequency g modes, mode coupling between g modes associated 
with different Xkm's might have significant influence on the frequency spectrum and the stability, as suggested by Lee (2001). 
For inertial modes in an isentropic star for whcih the Brunt- Vaisala frequency A'' vanishes identically in the interior, we found 
it difficult to obtain the oscillation frequencies in general. This is quite reasonable because the traditional approximation is 
considered to be valid for low frequency modes with N/uj >> 1. Interestingly enolugh, however, we do obtain the frequencies 
corresponding to the fundamental r modes with I' = m for isentropic polytropes even in the traditional approximation. The 
frequency ui of the r mode in the traditional approximtion is very close to 2n/(|m| + 1), which is consistent with previous 
calcualtions done without the traditional approximation (e.g., Yoshida & Lcc 2000a). 

Theoretically speaking, both g modes and r modes are excited in rapidly rotating SPB stars. As shown by Aerts et al 
(1999), a fraction of SPB stars are rapid rotators, and we think it is tempting to look for a link between Be stars and raidly 
rotating SPB stars. Recent discovery of a Be star that also shows SPB variability caused by g modes and r modes (private 
communication with G. Walker and with H. Saio, 2005) indicates such a possibility. If low m r modes excited in SPB stars 
accelerate the surface rotation effectively as a result of angular momentum redistribution in the interior, this acceleration 
can be a candidate mechanism for disc formation around the rapidly rotating star. We speculate unstable r modes are a 
steady supplier of angular momentum to the disc and prograde and retrograde g modes also excited in SPB stars may cause 
quasi-periodic pheonomena in Be stars. Further theoretical and observational studies are warranted. 
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Figure 1. Evolution tracks of main-scqiiciicc stars with mass ranging from 3Mq to 8Mq for X = 0.7 and Z = 0.02. The circles stand 
for SPB stars tabulated in Waelkens et al (1998), and the filled circles indicate rapidly rotating SPB stars as classified by Aerts et al 
(1999). 
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Figure 2. = co^./ \/ GMj B? versus LogTgff for odd r modes with = m = 1 of SMq main-sequence models for Q = Q/ ^JbmJb? 
0.4. The small thin circles and large thick circles stand for pulsationally stable and unstable modes, respectively. 
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Figure 3. Same as Figure 2 but for odd r modes with I' = m = 2. 
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Figure 4. Same as Figure 2 but for €1 = 0.2. 
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Figure 6. Growth rate rj = — oji/wr versus the inertial frame oscillation frequency \a\ of unstable odd r modes with I' = m = 1, 2, and 
3 for a 5M0 model with Log(Tefi) = 4.188, where n = 0.4 (0.0118mHz). The two filled wedges on the horizontal axis stand for Q and 
2n. 
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Figure 7. Same as Figure 2 but for even r modes with I' = m + \ = 2. No even r modes with I' = m + I = 2 are found pulsationally 
unstable. 
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Figure 8. Work integral w(r) of an even I' = ra+l = 2 r2o mode of a SM© main-sequence model with Log(Teff) = 4.188, where = 0.4. 



R modes of slowly pulsating B stars 15 



0.55 - 



0.45 - 




0.35 



4.26 



Log (T 



Figure 9. ojr versus LogTeft for odd r modes with I' = m = 1 oi 8Mq main-sequence models for Q = 0.6. The small thin circles and 
large thick circles stand for pulsationally stable and unstable modes, respectively. 
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Figure 10. Eigcnfunctions a;Rc(iT;/ ) versus the fractional radius x = r/R for odd I' = m = I r modes (panel a) and for even 

l' = m+ l = 2r modes (panel b) of the 5Mq ZAMS model, where fi = 0.4. In panel (a), the long dash-dotted line, short dash-dotted 
line, dotted line, and solid line represent tq, ri, rs, and rio modes, respectively. In panel (b), the long dash-dotted line, short dash-dotted 
line, dotted line, and solid line stand for ri, r2, re, and m modes, respectively. The amplitude normalization is given by Re(5(j) = 1 at 
the surface of the star. 
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Figure 11. Eigenfunctions xRe(iT;/^) versus x = r/R for odd V = m = \ r modes of a 5M0 model with Log(Tofi) = 4.188, where 

n = 0.4. Here, the dash-dotted line, dashed line, dotted line, and solid line represent ro, r\, rio, and r2o modes, respectively. The inlet 
shows ro mode. The amplitude normalization is given by Re(5(j) = 1 at the surface of the star. 
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Figure 12. t- ^{r) versus r/R for an odd V = m = 1 r20 mode of a 5Mq main-sequence model with Log(reff) = 4.188 for Q. = 0.4, 

where the solid line, dotted line, and dashed line stand for Tq^, t^^ , and r^^, respectively. The amplitude normalization to calculate 
the r20 mode is given by Re(5jj ) = 1 at the surface of the star. 
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Figure 13. T^^iK) and j versus a>R, for odd r modes with (' = m = 1 of a SMq main-sequence model with Log(Tg(5f) = 4.188 for 
?l = 0.4, where the circles and squares stand for Tq^{R) and J, respectively, and the open and solid symbols indicate pulsationally stable 
and unstable modes, respectively. The amplitude normalization to calculate the r modes is given by Re(5jj ) = 1 at the surface of the 
star. 
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Figure 14. iDr and growth rate \r}\ versus the radial order n oi I' = m = 1 r modes of the SMq ZAMS model for Q = 0.4. In panel (a), 
the small filled circles and large open circles stand for results obtained with and without the traditional approximation, respectively. In 
panel (b), the squares and circles stand for results obtained with and without the traditional approximation, respectively, and the filled 
and open symbols indicate pulsationally unstable and stable modes, respectively. 
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Figure 15. Same as Figure 14a but for even r modes with I' = m + 1 = 2. 
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Figure 16. Separation factor Xkrn versus the radial order ra of r modes with (k,m) = (—1, 1) (filled circles) and (—2, 1) (open circles), 
computed in the traditional approximation for the 5Mq ZAMS model for O = 0.4. The r modes associated with {k,rn) = (—1, 1) and 
(—2, 1) correspond to those with I' = m = 1 and I' = m + 1 = 2, respectively. 
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Figure 17. Eigenfunctions a;Re(5(j^) (solid line) and xRe(zi) (dashed line) are given versus x = r/R for the fundamental ro modes with 
I' = m = I (panel a) and for rio mode with I' = m = 1 (panel b) of the 5Mq ZAMS model, where Cl = 0.4. The eigenfunctions Si^ and 
zi arc computed without and with the traditional approximation, respectively. The amplitude normalization is given by Re(Sij) = 1 
and by Re(2;i) = 1 at the surface of the star. 
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Figure 18. Same as Figure 9 but for r modes computed in the traditional approximation. 



